Influence of the confinement geometry on surface superconductivity 
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The nucleation field for surface superconductivity, de- 
pends on the geometrical shape of the mesoscopic supercon- 
ducting sample and is substantially enhanced with decreasing 
sample size. As an example we studied circular, square, trian- 
gular and wedge shaped disks. For the wedge the nucleation 
field diverges as H c s/H C 2 = VS/a with decreasing angle (a) 
of the wedge, where H C 2 is the bulk upper critical field. 

PACS numbers: 74.80-g, 74.20De, 75.25Dw 



Recent progress in microfabrication techniques has 
made it possible to study mesoscopic superconducting 
samples of micronmeter and submicronmeter dimensions 
{{],§). The size and shape of such samples strongly in- 
fluences the superconducting properties. Whereas bulk 
superconductivity exists at low magnetic field (either 
the Meissner state for H < H c \ in type-I and type- 
II superconductors or the Abrikosov vortex state for 
H c i < H < H c2 in type-II superconductors), surface su- 
perconductivity survives in infinitely large bounded sam- 
ples up to the third critical field H c % sa 1.695£f C 2 0- 

For mesoscopic samples of different shape, whose char- 
acteristic sizes are comparable to the coherence length 
£, recent experiments have demonstrated the 

Little-Parks-type |Q] oscillatory behavior of the phase 
boundary between the normal and superconducting state 
in the H — T plane, where T and H are the critical 
temperature and the applied magnetic field, respectively. 
While for a circular sample, when the superconductor 
state is characterized by a definite angular momentum, 
the ocsillatory behavior of the nucleation field is well un- 
derstood theoretically p|-pd|, the problem of the nucle- 
ation of superconductivity in arbitrary shaped systems is 
still to be solved. 

Earlier numerical calculations by the present authors 
fllf showed that the oscillations of the nucleation field 
with sample size for ellipsoidal and rectangular shaped 
samples disappears only for large aspect ratios. Recently, 
the asymptotic behavior of the nucleation field in the 
small size limit, L £, and in the large size limit, L -C £, 
where L is the characteristic sample size, was found for 
smooth (i.e. without sharp corners) and square samples 
| p2[ . The particular case of surface superconductivity in a 
rectangular loop [fl3f and in a wedge [ 14| was investigated 
recently. In Ref. p4| a variational approach was used and 
they found that H C 3 is maximum for a wedge angle of 
a = 0.447T. This is surprising in view of the intuitive idea 
that surface superconductivity increases with decreasing 



size of the system and therefore, one would expect that 
H C 3 is an uniform increasing function with decreasing 
a. The result of Ref. jlij was also in disagreement with 
calculations [ pd|Jl2| for the particular case of a square 
sample. 

In the present paper, we calculate, within the 
Ginzburg-Landau (GL) mean field approach, the nucle- 
ation field for bounded square and triangular samples as 
well as for a wedge with infinitely long sides and arbitrary 
corner angle a. We resolve the discrepancy between the 
results of Refs. jl3] and [E3. Furthermore, we obtain the 
following analytical expression for the asymptotic behav- 
ior of the nucleation field of a wedge H C 2,/H C 2 = y/3/a 
in the limiting case of a <C 1. 

Nucleation of superconductivity in a finite bounded 
sample. Close to the superconducting - normal transi- 
tion the demagnetization effects are not important and 
consequently the nucleation of superconductivity is de- 
scribed by the first GL equation for the order parameter 
<F 
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where a, (3 are the GL parameters which depend on tem- 
perature, A is the vector potential of the applied mag- 
netic field H = rot A and m* = 2m and e* = 2e are the 
mass and electrical charge of the Cooper pairs, respec- 
tively. We consider a uniform external magnetic field 
applied along the z— axis and assume that the order pa- 
rameter is uniform in this direction. The latter is sat- 
isfied for samples which are close to the nucleation field 
and which have flat top and bottom sides. When calcu- 
lating the nucleation field we neglect the nonlinear term 
and rewrite Eq. (m) as 
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where the distance is measured in units of the coherence 
length £ = h/y/— 2ma, the vector potential in cS./2e£, 
and the magnetic field in H C 2 — ch/2e^ 2 , the bulk upper 
critical field. 

Thus the problem of nucleation of superconductivity is 
reduced to the calculation of the lowest eigenvalue and 
corresponding eigenfunction of the two-dimensional (2D) 
linear operator L = (— iV — A) 2 which also describes 
the motion of a quantum particle. The transition from 
the normal to the superconducting state occurs when the 
lowest eigenvalue of L becomes smaller than unity. Con- 
trary to the case of the usual Schrodinger equation, where 



1 



the wave function is taken equal to zero at the sample 
boundaries, here the normal component of the supercon- 
ducting current must be zero at the boundary between 
the superconducting and the insulating material 



iV- A) Ltf = 0. 



(3) 



This difference in boundary condition changes drastically 
the behavior of the eigenvalues as function of the sam- 
ple size. Whereas the energy of a quantum particle in- 
creases as E ~ 1/L 2 for decreasing sample size L, the 
resulting energy for the boundary condition Eq. (3) de- 
creases as E ~ L 2 . As a result the nucleation field de- 
creases with increasing sample size and tends to the limit 
_ff C 3 = 1.695i?c2 for infinitely large samples with a flat 
boundary. In circular samples this increase of H c % with 
decreasing sample size has an oscillatory behavior due 
to the quantization of the angular momentum. For ax- 
ially symmetric samples the latter is an exact quantum 
number 

To find the nucleation field for an arbitrary shaped 
sample we use a finite-difference representation of the 
differential operator L on a uniform Cartesian space grid 
within the link variable approach Jl5[ . The lowest eigen- 
value and corresponding order parameter are found by 
applying the following iteration procedure L^ 1 = 
The calculated nucleation field as a function of the square 
root of the sample area is shown in Fig. 1 for different 
geometries. Note, that the oscillatory behavior of the 
nucleation field survives even in triangular samples al- 
though the oscillations are less pronounced then for a cir- 
cular geometry. The contour plots of the order parameter 
density |^| 2 are depicted in the insets of Fig. 1 for two 
square samples with different sizes. For smaller samples 
areas S — (2.32£) 2 , the maximum of the order parame- 
ter is situated in the center of the square, which corre- 
sponds to the Meissner state (see upper inset of Fig. 1). 
The jump in the first derivative of the nucleation field at 
y/S w 2.33£ is caused by the appearance of a vortex in 
the center of the square (see lower inset of Fig. 1). This 
changes the slope of the nucleation field and corresponds 
to the appearance of a vortex state with a larger value 
of the angular momentum. The existence of a sequence 
of vortex-like states is inherent to all considered samples. 
Notice also that the nucleation field is larger for samples 
with sharper corners as expected intuitively. 

The nucleation of superconductivity in a wedge. For 
the wedge geometry it is more convenient to use cylindri- 
cal coordinates (p,<f>) with = Hp 2 /2 and to measure 
all distances in */ ch/eH . Then the linearized first GL 
equation takes the form Jl4[ : 
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The nucleation field H c ^ — 2H C 2/X is obtained from the 
lowest eigenvalue A of the operator in the LHS of Eq. (Q). 
In order to find the lowest eigenvalue and the correspond- 
ing eigenfunction we use again the finite-difference tech- 
nique described above. Instead of an infinitely long wedge 
we consider a finite fragment with a very large radius 
R = 15/y/a such that the value for A is independent of R. 
Fig. 2 shows the spatial distribution of the square of the 
absolute value of the order parameter \^f\ 2 for two differ- 
ent wedge angles a = 0.057T and 0.57T. For the a = 0.05n 
case, the order parameter practically does not depend on 
the azimuthal angle and decays faster then exponentially 
deep into the sample. For a = n/2, the order parameter 
still decays quickly with the radius but a prominent an- 
gular dependence appears, specially for large radii. This 
expected behavior, namely superconductivity only exists 
in the corner [fl3| , implies a decay of the order parameter 
from the wedge side into the sample. 

The numerically obtained nucleation field is shown in 
Fig. 3 by the solid curve which decreases monotonically 
with increasing wedge angle and diverges in the limit 
a — ► 0. For a = ir/2, we found H C 3 1.96i/ C 2 which 
is close to the estimated result H c ^ w \.S2H C 2 of Jl2| . 
Note, that increasing the wedge angle beyond a > it/2 
changes very weakly the nucleation field and we found 
the well-known result H c ^ — 1.695_ff C 2 for a — it |l6| . 

Because the order parameter varies only weakly with 
the angle it allows us to find the asymptotic behavior of 
the order parameter analytically for small wedge angles 
a -C 1. To this end we rewrite Eqs. (||J|) as 
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where x = \fa.p, i] — <p/a, p — A/a, ip — exp{—ix 2 rj)^! 
Ip^l . The second term in the LSH of Eq. (|^) dominates for 
the case of very small wedge angles. Therefore, the new 
order parameter ip depends slightly on the angle rj which 
is also in agreement with the boundary condition (0). 
This allows to simplify the problem. Using the boundary 
condition (Q) we integrate Eq. (^) over ry assuming that 
ip depends only on the radius and obtain 



d 2 ip 1 
dx 2 x 



. c^ip .„ Ax 2 . , 
2^-(2*--)<A 



pip. (8) 



Note, that the wedge angle does not enter into the last 
equation and the reduced eigenvalue p does not depend 
on a. After the substitution ijj(x) = exp(—ix 2 /2)f(x) 
Eq. (||) is reduced to the well-known equation for the 
harmonic oscillator with the lowest eigenvalue p = 2/ 
Thus the order parameter for sharp wedges a <C 1 can 
be written as 
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¥ = aBp((*0-i| ~^)P 2 ) . ( 9 ) 

and the nucleation field is inversly proportional to the 
wedge angle 

H c3 = ^-H c2 . (10) 
a 

Our numerical results deviates from the asymptotic ex- 
pression ( |io| ) with about 10% when the wedge angle is 
increased up to a ~ 0.15-7T w 0.5 (compare the dashed 
curve with the solid curve in Fig. 3). As is evident from 
Eq. (^|), the corrections to the above asymptotic nucle- 
ation field are of second order 0(a 2 ). From this obser- 
vation we tried to fit H c3 /H C 2 to a function g(a 2 )y/3/a. 
Within the accuracy of our numerical results we found 
that the nucleation field for a wedge could be accurately 
fitted to 

H c3 V3 / , 0.746a 2 \ , x 

77^ = — 1 + 0.14804a 2 + ^ — — . 11 

H c2 a V « + 1.8794/ 

This function is shown by the dotted curve in Fig. 3 and 
agrees very well with our numerical results (solid curve). 

Let us compare the asymptotic result for the nucleation 
field in a wedge with those for thin strips and small cir- 
cles. Assuming that the order parameter depends slightly 
on the width d of the strip and on the radius R of the 
circle, we find the following asymptotic expressions for 
the nucleation field in the strip H c3 = V3H C 2£,/d and the 
circle H c s = 2\/2H C 2^/R. In conlusion, for L -C £, where 
L is the smallest dimension of the sample, the nucleation 
field increases inversely with the sample size, which cor- 
responds to a£ in the case of a wedge-like sample. 

Note added in proof: The authors of Ref. G3] re- 
cently found a mistake in their calculation which was cor- 
rected in the erratum |17|| . The new asymptotic behav- 
ior found with their variational calculation is H c % / H C 2 = 
= 0.816/a which is a factor 3/V2 = 2.12 smaller 
than our exact result (10). 
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FIG. 1. The nucleation field as function of the square 
root of the sample area S for different sample geometries. 
The insets show |\t| 2 for square samples with different side 
lengths \f§. 

FIG. 2. The distribution of \&| 2 in a wedge sample for 
two different wedge angles a = 0.057T (a) and a = 0.57T (b). 

FIG. 3. The nucleation field in a wedge sample as function 
of the wedge angle a. Solid curve shows the results of our 
finite-difference calculation, the dashed curve corresponds to 
the asymptotic result for small angles and the dotted curve is 
the fitted Eq. (11). The open and solid circles are results for 
a = 7r/2 [12] and a — n/2 [16], respectively. 
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